On the structure of blue phase III 
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We report large scale simulations of the blue phases of cholesteric liquid crystals. Our results suggest a 
structure for blue phase III, the blue fog, which has been the subject of a long debate in liquid crystal physics. 
We propose that blue phase III is an amorphous network of disclination lines, which is thermodynamically and 
kinetically stabilised over crystalline blue phases at intermediate chiralities. This amorphous network becomes 
ordered under an applied electric field, as seen in experiments. 
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Liquid crystals (LCs) offer prime examples of systems with 
spontaneously broken symmetry that support topological de- 
fects of various types Q. Alongside their technological util- 
ity in various optical devices, LCs offer testing grounds for 
fundamental theories whose counterparts in other fields (such 
as cosmology, particle physics, and more recently exotic su- 
perconductivity Q) are less directly testable. While many 
aspects of LC physics are by now understood, some, such as 
the character of the so-called 'blue fog' (blue phase III) have 
remained unresolved despite efforts spanning several decades. 

The simplest liquid crystalline phase is the nematic, in 
which molecules have a preferred orientational axis (the di- 
rector) but no translational order. Introduction of molecular 
chirality causes the director to precess in space. In the sim- 
plest resulting phase (the cholesteric) it does so about a single 
helical axis, thereby creating a 1-D periodic structure whose 
wavelength is the helical pitch. Locally, however, the ordering 
remains nematic and the cholesteric therefore supports topo- 
logical line defects, known as disclinations, in which the di- 
rector executes a rotation of tt under a full 2tt rotation around 
the defect line. Within the cholesteric phase itself, these exci- 
tations are absent from the ground state. 

In many chiral liquid crystals the transition between 
the cholesteric and isotropic phases occurs through a cas- 
cade of weakly first-order transitions to intermediate struc- 
tures, known as 'blue phases' (BPs) and consisting of self- 
assembled disclination networks (BE)- At high enough chi- 
rality, a simple helical structure is less stable locally than a 
so-called double-twist cylinder (DTC), in which the director 
field rotates simultaneously in two directions perpendicular to 
a cylinder axis. However DTCs cannot be smoothly patched 
together to fill the whole of space; disclination lines are re- 
quired at the interstices between the cylinders. If an external 
field is used to align the DTCs, the disclinations can form a 
simple linear array |4], but more generally they meet at junc- 
tions forming a multiply connected network. The resulting 
structure is highly colored ('blue') since the spacing between 
defects lies in the optical range. BPs are stable when the free 
energy gained by creating DTCs instead of a simple helix is 
enough to compensate the free energy loss due to defect for- 
mation. They therefore occur at high chirality (where DTCs 
are favored) and for T close to but below Tic (the transition 
temperature to the isotropic phase) where defect energies are 
low. Long viewed by many as little more than a curiosity Q, 



BPs have recently emerged as promising materials for pho- 
tonics and display devices following their stabilisation over a 
much larger temperature range than previously thought pos- 
sible (about 60 K compared to <1 K) 0. This has revived 
interest in these fascinating hierarchical materials, in which 
molecular-scale interactions stabilize topological defects that 
self-organize into a micron- scale superstructure. 

Three BPs have been observed experimentally (at zero elec- 
tric field). Two, BPI and BPII, are highly ordered: their discli- 
nation networks form a regular cubic lattice, and their physics 
is well understood. (Another ordered phase, O5 , is predicted 
by theory at high chirality, but not found experimentally.) The 
last one is known as BPIII, and its structure is one of the un- 
solved puzzles of liquid crystal physics. Theorists have pro- 
posed that BPIII may either be a quasicrystal 0, a spaghetti- 
like tangle of double-twist cylinders |7 |, an amorphous state 
formed by BPII domains 0, or a metastable phase |9|. Ex- 
periments on the thermodynamics, scattering and electric field 
response [ T0HI2I remain inconclusive, although electron mi- 
crographs appear to favour an amorphous structure and to rule 
out the quasicrystal model 0. 

In this work we show that computer simulations can help 
settle this important physical question. By means of large 
scale (supra-unit-cell) simulations, which enable very accu- 
rate comparison of free energies for both ordered and amor- 
phous structures, we provide strong evidence that BPIII is in- 
deed an amorphous network of disclinations. Remarkably we 
show that, within a certain window of chirality and with a 
standard choice of free energy functional detailed below, in- 
dividual aperiodic structures exist that are more stable than 
either BPI, BPII, or O5 . This narrow window lies within the 
observed range of BPIII, and can only become wider, in line 
with experiments, if one allows for configurational entropy. 
(Such entropy arises when there exist multiple aperiodic struc- 
tures of the same free energy, as seems likely here.) Further- 
more, we show that an applied electric field orders our aperi- 
odic BPIII candidate into a different, much more ordered blue 
phase. This concurs with longstanding experiments [Qj] |T2) 
which showed evidence of a field-induced transition to a new 
phase, BPE, whose structure was not previously identified. 

The local order in the BPs can be described by a trace- 
less, symmetric, second rank tensor, Q whose equilibrium 
thermodynamics is governed by the Landau-de Gennes free 
energy functional T, whose form, within the one-elastic con- 
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stant approximation (discussed later) is standard, and speci- 
fied at fT3l . Within that approximation, the phase behavior 
of BPs depends on just three dimensioness parameters, a re- 
duced temperature r, a reduced chirality ft and a reduced field 
strength £ lfT4l . Expressions for these in terms of K (the 
elastic constant), qo (the cholesteric wave- vector) and other 
parameters in T are given at fT3ll . Good agreement between 
theoretical 1 15 ] and experimental 1 16 ] phase boundaries is ob- 
tained by taking r ex (T — Tjc) and ft linear in the mole- 
fraction of a chiral component (with respective proportionality 
constants ~ 2 K _1 and ~ 2 for one specific mixture (El). 

We employ a 3D hybrid lattice Boltzmann (LB) algo- 
rithm fT5l [T7UT91 to solve the Beris-Edwards equations. The 
evolution of the Q tensor l20l is 

AQ = r(^f + i Tr (|§)i). (i) 

Here, T is a collective rotational diffusion constant and D t is 
a material derivative for rod- like molecules |2Q| . The term 
in brackets is the molecular field, H, which ensures that Q 
evolves towards a minimum of the free energy. The fluid ve- 
locity field obeys the continuity equation and a Navier-Stokes 
equation with a stress tensor generalised to describe liquid 
crystal hydrodynamics, and discussed elsewhere ifTTl . 

Though these equations represent the true dynamics, we 
use them here simply to find free energy minima. Thus, as 
in previous work we additionally allow a so-called 'redshift' 
in which the parameters in T are dynamically updated at fixed 
r, ft, £ fT3l . This exploits a scaling among those parameters to 
ensure that the system is not frustrated by periodic boundaries: 
in particular, for any cubic BP, a lattice parameter emerges 
that truly minimizes T. The accuracies of their computed free 
energies F = min(J r ) are (at least for £ = 0) thus limited 
only by discretization. This is chosen to fully resolve the de- 
fects fT3l . whose core energy is finite, and set by T itself. 

In computing the free energy of periodic structures (BPI, 
BPII and O5) we apply a perturbation of the appropriate sym- 
metry to a uniform state and then evolve dynamically fT3l . 
This delivers an accurate F value in each such phase, of which 
the lowest can be chosen; but (in common with all other meth- 
ods for computing free energies of ordered phases) we cannot 
rule out others of still lower T. To address BPIII, we also 
need to generate aperiodic candidates. Here our strategy is 
similar: we start from various different aperiodic initial con- 
ditions, evolve each dynamically, and choose that of lowest 
T |[T3ll . As shown below, this beats all three periodic struc- 
tures within a certain parameter window. Because we cannot 
exhaust all possible initial conditions, our free energy is an 
upper bound on aperiodic states; further exploration can thus 
only widen that stability window. In practice, all the aperi- 
odic candidates we generated look similar: fuller exploration 
is thus unlikely to change our conclusions about the character 
of BPIII. 

To minimize finite size effects, we simulate very large sys- 
tems (in contrast to (HJ). Typically we used 128 3 lattices, 
which accomodates 8 half pitches in each direction l2T1l . Se- 
lected simulations with 25 6 3 lattices confirm the results which 
we report below. As shown in fT9lL it is easy to generate ape- 
riodic structures by placing a localized nucleus of BPI or II 



in a cholesteric or isotropic matrix. However, the lowest T 
values we have so far found are instead achieved by initial- 
ising the system in the cholesteric phase in the presence of a 
low density (typically about 1-2 % in volume) of randomly 
placed doubly twisted droplets. Once initialized, the system 
is relaxed dynamically until it reaches a quiescent end- state. 

For low chirality (ft < 1.5), the initial defects are washed 
out to leave a cholesteric phase. However, for a large regime 
of intermediate chiralities (1.75 < ft < 3), our simulations 
show an intriguing dynamics, through which the dilute dou- 
bly twisted regions grow and rearrange dramatically to form 
a whole network of disclinations, which very slowly creeps to 
an amorphous end- state (Fig. 1, top). Its amorphous charac- 
ter is confirmed by the structure factor C(k) (Fig. 1, bottom). 
The ring in C(k) is set by the average distance between the 
branch-points in the defect network, whereas the small resid- 
ual peaks (which break spherical symmetry) are likely due to 
residual finite size or periodic boundary effects. The amor- 
phous network approaches equilibrium through very slow lo- 
cal rearrangements of the disclination junctions; the end-state 
in Fig. 1 is very close to kinetic arrest, in at least a local mini- 
mum of T. Notably, four disclination lines meet at most junc- 
tions, so that, as suggested by experiment 1 10 ], the structure 
is locally closer to that of BPII than BPI. 

Remarkably, for larger values of the chirality (ft > 3) our 
simulations attain a much more regular state, closely resem- 
bling 5 which (see below) minimizes our chosen T at very 
high ft. (In O5 itself, whose free energy we have computed 
precisely, eight disclination lines merge at each junction fT3l .) 
Therefore, our methodology is capable in principle of finding 
a periodic disclination lattice, if kinetically accessible. Given 
that our equations fully incorporate liquid crystal hydrody- 
namics, we believe that the kinetic propensity to form a disor- 
dered disclination network for intermediate chiralities indeed 
reflects a real physical property, although our limited simu- 
lation times (~ 1ms fT3ll ) may well exaggerate the stability 
window of amorphous structures in the (ft, r) plane. 

A crucial question is whether our BPIII candidate structure 
is only kinetically, or also thermodynamically stable. To an- 
swer this, we have carefully compared F(k,t) for BPIII with 
those of BPII and O5, for a range of chiralities ft at selected 
values of r. (BPI is not competitive in the ft range of interest 
here.) As an example, curves for F(k, —0.25) are shown in 
Fig. 2. We see that there is a small but finite chirality window 
in which the BPIII-network is the thermodynamic equilibrium 
phase. We find it remarkable that any single aperiodic struc- 
ture can outcompete such periodic ones when minimizing the 
relatively simple Landau-de Gennes free energy. In partic- 
ular, this minimization takes no account of order-parameter 
fluctuations about the local minimum, which might help sta- 
bilize BPIII l22l . nor the configurational entropy associated 
with having many such minima. (These neglected contribu- 
tions should be small - of order ksT per unit cell while the 
free energy differences in Fig. 2 are typically of order 100- 
1000 times larger fT3l .) However, the thermodynamic stabil- 
ity window in BPIII is smaller than the one we find kinetically 
- in agreement with the expectation that amorphous struc- 
tures should form more easily dynamically than highly or- 
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FIG. 2: Free energy densities (V is the volume) vs. chirality k at r — 
—0.25 for BPII, our candidate for BPIII, and O5. The inset shows a 
blow-up of the region where our BPIII is the equilibrium phase. Our 
free energy densities are expressed in units of Kq$ (see 1 13 1 for the 
mapping between simulation and real units). 



FIG. 1: Top: End-state disclination network at r — —0.25, k — 2.5: 
The picture shows the isosurface q(r) = 0.12, with q the largest 
eigenvalue of Q. Inside each tube is a disclination line (on which q 
takes a minimum value). Bottom: Structure factor C(k) = |g(k)| 2 , 
on cuts along k y = (left) and k x — (right) with wavevectors 
k x /qo,k z /q £ [-4, 4] and k y /q ,k z /qo G [-4, 4] respectively. 



dered counterparts, which require long and complicated pro- 
cess to annihilate any dislocations and overcome relatively 
high energy barriers. The window of thermodynamic stability 
might also be broadened by relaxing the one elastic constant 
approximation, particularly if this raises the free energy of O5 
relative to the other states. This certainly merits further study, 
given the experimental absence of that phase. 

All the above arguments support our assignment of the 
amorphous network seen in Fig. 1 as the theoretically elusive 
BPIII, stable at larger k than BPII, but locally similar to it 0. 
Experiments on BPIII in electric fields also point to a field- 
induced transition between BPIII and another structure which 
has been named BPE HUEl. BPE was found to give rise to a 
sharper peak in the scattering data consistent with an enhanced 
ordering (unfortunately Refs. |[TT1 [TZ l do not provide further 
experimental data on the structure of BPE). It is therefore in- 
teresting to ask what happens when our amorphous structure 
is subjected to an electric field. In our simulations we can 
follow the evolution of the Q tensor and of the disclination 
network in an external field £ , and also compute C(k); if a 
crystalline phase emerges, it will exhibit Bragg peaks. 

By stepwise increasing £ from to 0.65, we found that the 
disclination network in Fig. 1 melts away leaving a nematic 
state. (This holds for positive dielectric anisotropy [13].) We 
then performed a long simulation at £ = 0.55, close to but be- 
low the threshold beyond which the system becomes nematic. 



The time evolution of our BPIII, stable at zero field, is shown 
in Fig. 3. This shows an interesting rearrangement of the de- 
fect texture, through which the junction points rotate and fi- 
nally reconstruct to yield a topologically distinct phase, with 
helical disclinations lying along layers stacked perpendicular 
to the electric field (which is vertical). The disclination lines 
in two consecutive layers are turned by 90 degrees, so that 
they show a square arrangement when viewed along the field 
direction. Within a layer two adjacent disclination lines are 
staggered. Again in agreement with experiments, then, we 
find that our amorphous BPIII candidate undergoes a field- 
induced transition to an ordered disclination network which 
we provisionally identify as BPE. It would be interesting to 
perform additional experiments to test such an identification. 
Our candidate BPE appears related to, but distinct from, other 
field- or confinement-induced BPs previously reported l4l l23lL 
with strong crystalline order (albeit with some defects remain- 
ing) as confirmed by its structure factor. 

In conclusion, we have presented large scale simulations of 
the cholesteric blue phases. We have given strong evidence 
that in a region of the temperature-chirality plane, instead of 
a periodic BP, an aperiodic one is selected kinetically. More- 
over, and strikingly, within the (one-elastic constant) Landau - 
de Gennes free energy fl3lL we have found a finite window of 
chirality in which this amorphous network is thermodynami- 
cally more stable than the competing crystalline blue phases, 
BPI, BPII and O5. These facts suggest that our kinetically 
and thermodynamically stable network is none other than the 
blue fog, BPIII. This view is strengthened by the study of the 
field response of this structure, which reconstructs into an or- 
dered phase at intermediate values of the electric field, slightly 
smaller than those at which the disclination network melts to 
give a field-oriented nematic phase. 

We believe our simulations shed important light on the 
structure of the blue fog, which we propose is an amorphous 
disclination network which is locally close to BPII: this view 
is perhaps intermediate between the "BPII domain" and 
the "double twist tangle" [ 7 ] models which were proposed in 
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FIG. 3: Isosurfaces of q = 0.12 at r = 0.2, k = 2.0, S = 0.55: The 
sequence shows the transition from the amorphous BPIII network to 
a field-induced BPE. The field direction, z, is vertical. Times shown 
are t ~ 0.15,0.3,0.6 and 1.4 ms (taking one timestep ~ 1 ns). 
Bottow row: Structure factor C(k) at the end of the run. The pictures 
show a cut through the 3D data at k y = (in the (k x ,k z ) plane, left) 
and k z = (in the (k x , k y ) plane, right), with scales as in Fig. 1. 



the early 1980's, and is at odds with some of the subsequent 
structures proposed for this phase. Several open questions re- 
main or are stimulated by our results. Firstly, one should char- 
acterise statistically the phase transition from the blue phase 
to the isotropic phase: this may be possible by adding ther- 
mal noise to simulations like the ones reported here. Sec- 
ondly, use of more accurate scattering and visualisation tech- 
niques should allow tests of our candidate structure for BPE 
(see Fig. 3). Finally, if the large energy scales that separate 
the various BP topologies (of order 100 — 1000/c^T per unit 
cell fT3ll ) also control reconstruction of the defect network 
within BPIII itself, one may expect that in some materials, 
even when thermal noise is allowed for, such reconstruction 
cannot be achieved on any reasonable timescale. If so, BPIII 
will represent another elusive entity: an amorphous state of 
globally minimum free energy that is kinetically arrested — 
or an 'equilibrium glass'. 

We thank A. N. Morozov for useful discussions. This 
work was supported by EPSRC grants EP/E045316/1 and 
EP/E030173/1. MEC is funded by the Royal Society. 



Here we provide further information on the free energy 
functional and the control parameters r, k,£ (section A); on 
resolution and on 'redshift', i.e., the parameter-dependent unit 
cell size of blue phases (section B); on the initial conditions 
chosen for runs presented in the main text (sections C and D); 
on the conversion of simulation parameters to physical units 
(section E); on the free energy comparison between different 
blue phases (section F); and on the effect on field-response of 
changing the sign of the dielectric anisotropy (section G). 



1. Free energy functional, chirality and reduced temperature 

The thermodynamics of cholesteric blue phases can be 
described via a Landau-de Gennes free energy functional 
T[Q] = / <i 3 rf(Q(r)). The free energy density we use, fol- 
lowing [ 1 ], is: 

/(Q) = ^(l-l)Ql -^f Qa ,Q 01 Q la+ ^L { Ql^ 
K K 

- j^E a Q a pE(3 (2) 

Here repeated indices are summed over and Q^a stands for 
QapQap etc. The first three terms are a bulk free energy den- 
sity whose overall scale is set by Aq (discussed further be- 
low); 7 is a control parameter, related to reduced temperature. 
Varying the latter in the absence of chiral terms (qo = 0) gives 
an isotropic-nematic transition at 7 = 2.7 with a mean-field 
spinodal instability at 7 = 3. The next two terms in Eq|2] 
describe distortions of the order parameter field. As is con- 
ventional d [24| we have assumed that splay, bend and twist 
deformations of the director are equally costly; K is then the 
one elastic constant that remains. The parameter q Q is related 
via go = 2n/po to the pitch length, po, describing one full 
turn of the director in the cholesteric phase. 

The remaining term accounts for the coupling to the electric 
field, E a , and e a is the dielectric anisotropy. In the text we 
considered the case s a > 0, whereas below we show what 
happens for negative dielectric constant. 

The phase behaviour of blue phases (in the absence of noise 
due to thermal fluctuations, i.e., as computed by minimization 
of J 7 ) may be shown to only depend on the following dimen- 
sionless parameters Q : 

r = 27(1 - 7 /3)/ 7 (3) 
K = yj 108 Kql/(A Ql ) 
£ 2 = (27s a E a E a )/(327rA 0l ). 

as used in the main text. These are referred to as reduced tem- 
perature, chirality and reduced field respectively. Note that if 
/ is made dimensionless, r appears as pref actor of the term 
quadratic in Q, whereas k quantifies the ratio between bulk 
and gradient terms. 
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2. Resolution, unit cell size and 'redshift' 

In our simulations, each disclination core spans a few (typ- 
ically 3-4) lattice sites. To confirm that the defect core is ade- 
quately resolved, we have run selected simulations in which 
the resolution was doubled or quadrupled, and found very 
similar values for the free energy. On this basis we assess 
that the typical error due to discretization for our main runs is 
much smaller than the free energy difference between BPs. 

In cubic BPs it is known that the equilibrium unit cell size is 
slightly different from (and typically slightly larger than) the 
half pitch of the cholesteric liquid crystal, ir/qo. To account 
for this, a 'redshift' factor r is introduced 031 whose variation 
effectively allows free adjustment of the BP lattice parameter, 
A — > A/r, despite the use of periodic boundary conditions in 
our simulations. That is, to avoid changing the size of the sim- 
ulation box, redshifting is performed by an equivalent rescal- 
ing of the pitch parameter and elastic constant, qo — » qo/r 
and K — >• K r 2 . In simulations aimed solely at free energy 
minimization, it is legitimate to make r a dynamic parameter 
and update it on the fly to achieve this. At each time step one 
then selects the value of r for which the physical free energy 
density is lowest. More details on this procedure are given in 
Ref. 0311 where this scheme was first proposed. 

We have used this dynamic updating of the redshift (i) in all 
single unit cell runs used to map out the free energy curves, 
and (ii) in multi-unit cell dynamical runs after the final net- 
work was formed. Redshifting is essential in order to get ac- 
curate estimates for the free energy of the cubic blue phases, 
whose relative free energy differences between phases is less 
than a percent of their typical values. With redshifting, the 
accuracy attainable is in principle limited only by discretiza- 
tion error; without it, the dominant error is from mismatch 
between the preferred unit cell size and the periodic bound- 
ary conditions. For amorphous networks, those conditions en- 
force a periodic structure (on the box scale) with the aperiodic 
structure within that; redshifting likewise avoids any incom- 
mensurability between the two length scales. Because of the 
constrained periodicity, the computed free energy for amor- 
phous structures will always be overestimated, despite which 
we find a window of thermodynamic stability for these with 
free energies lower than the known ordered phases. (As usual 
we cannot rule out a still lower free energy for some other, un- 
known, ordered phase.) Our use of redshift on all phases thus 
means that the true thermodynamic stability window of BPIII 
can only be wider than the one found here. 



3. Initial conditions: cubic blue phases 

To accurately compute the free energy of cubic blue phases 
from single unit cell simulations, we must start from an initial 
condition of the correct symmetry and disclination line topol- 
ogy. To this end, as discussed in Ref. lfT4l a. we can adopt 
the form for Q(x, y, z) corresponding to the limiting form for 
each phase in the high chirality limit d. 



Qxx — 
Qxy 
Qxz 



Q, 



yy 



For BPI, with q' = y/2q , we then have initially: 

" ~ ' -2 cos(%y) s'm(%z) + sin (%x) cos (%z) 
+ cos (q f x) sin (q f y) 

cos (q' y) cos (q' z) + y/2 sin (q f x) sin (q' Q z 

- sin (q' x) cos (q f y) 

~ \/2cos(q' x) cos(q f y) + y/2 s'm(q' z) sm(q' Q y) 

- cos(%x) sm(q z) 
—2 sm(^Qx) cos(q' z) + sm(q f y) cos(q f x) 

+ cos(q f y) sin(q' z) 

~ y/2 cos(q z) cos(q x) + y/2 sm(q f y) sin(q x) 

- sin((7o?/)cos((7o^). 

Likewise, for BPII we start from: 

Q xx ~ cos(2q z) - cos(2q y) 

Q xy ~ sm(2q z) 

Q xz ~ sin(2q y) 

Q yy ^ cos(2go^) — cos(2go^) 

Q yz ~ s'm(2q x). 



yz 



(4) 



(5) 



Finally, for 5 we start from: 



Q xx ~ 2cos(y/2q y)cos(y/2q z) 

— cos(y/2q x) cos(y/2q z) 

— cos(y/2qox) cos(y/2qoy) 
Q xy ~ y/2 cos(y/2q y) sm(y/2q z) 

— y/2cos(V2qox) s'm(y/2qoz) 

— sm(y/2q x) sin(y/2q y) 
Qxz - V2cos(V2q x) s'm(V2q y) 

— y/2 cos (V2 q$ z) sin(y/2 q^y) 

— sin(A/2<7o^) s'm(y/2qoz) 
Qyy ~ 2 cos(V / 2go^) cos(y/2qoz) 

— cos(y/2q y) cos(y/2q x) 

— cos(y/2q y) cos(y/2q z)) 
Q yz ~ V / 2cos(V / 2<7o^) sm{V2q^x) 

— y/2 cos(y/2qoy) s'm(y/2qox) 

— s'm(y/2qoy)s'm(y/2q z). 



(6) 



This procedure was already discussed in Refs. [15] and 
fT4l a, where additional details can be found. Fig. 4 shows the 
unit cell of O5 obtained after a full numerical minimisation 
(with redshift) starting from the high chirality limit detailed 
above. To ensure that equilibrium was reached, before ending 
a run we required that the free energy did not change, in two 
successive measurements, by more than an accuracy threshold 
(typically chosen as less than 0.1%). 

In line with previous theoretical and numerical calcula- 
tions dEL we find that O5 is the equilibrium structure at 
high enough chirality (see Fig. 2 in the main text). This is 
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FIG. 4: Image of the unit cell of O5 after free energy minimization. 



therefore expected within the free energy functional formal- 
ism that we use. Given that O5 is not observed experimentally, 
it would be interesting to study its window of thermodynamic 
stability by relaxing the one elastic constant approximation 
employed in this work. 



4. Initial conditions: blue phase III 

In the main text we presented configurations for our candi- 
date blue phase III which were obtained by initialising the sys- 
tem in the cholesteric phase in the presence of a low density 
(typically about 1-2 % in volume) of randomly placed doubly 
twisted droplets. One may therefore wonder what happens 
when different initial conditions are selected. 

To address this question, we have performed additional runs 
starting from a different initial configurations, in which the 
system was initialised in the isotropic phase (rather than in 
the cholesteric phase as in the main text), again in the presence 
of a low density of randomly placed doubly twisted droplets. 
Fig. 5 shows the initial and final disclination network in this 
case - it can be seen that the final state is very similar to the 
one shown in the main text. A minor difference is that, in the 
case of an isotropic background, we found that the initial vol- 
ume density of doubly twisted impurities needs to be slightly 
larger for the network to grow. In no case did we find a free 
energy by this route that was lower than the one resulting from 
the cholesteric initial condition. As shown in [19], finally, it 
is also possible to generate locally stable aperiodic structures 
by placing a localized nucleus of BPI or II in a cholesteric or 
isotropic matrix. These were found to have higher free ener- 
gies than those reported here (but comparable with them). 

To summarise, the kinetic stabilisation of amorphous net- 
works occurs from a wide range of initial conditions, includ- 
ing doubly twisted impurities inside a cholesteric (main text) 
or isotropic (background), and cubic blue phase nuclei inside 
cholesteric or isotropic phases. Such a predicted stabilisation 
is therefore a very robust result of our simulations. 



5. Parameter mapping to physical units 

Here we describe how simulation parameters are related 
to physical quantities in real BP materials. In order to get 



from simulation to physical units, we need to calibrate scales 
of length, energy and time. We follow the methodology 
of (17] [I9ll . First we define a set of lattice Boltzmann units 
(LBU) in which the lattice parameter £, the time step At, and 
a reference fluid mass density po are all set to unity. This 
is the set of units in which our algorithm is actually written. 
The length scale calibration is straightforward, and fixed by 
the cholesteric pitch po, which is typically in the 100-500 nm 
range 1 1 ]. More precisely, in our simulations we set the unit 
cell of BPI/II to be 16 LBU; this gives good resolution with- 
out wasting resource. Therefore the LBU of length (one lattice 
site) corresponds to, say, lOnm in physical space. 

To get an energy scale, we use the measurements cited in 
Appendix D of 1 1 ] which, as shown in |fT9l suggest: 
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2 - 5 x 10" b J 



(7) 



From this relation (given that 7^3 close to the 
isotropic/cholesteric transition) we obtain that Aq ~ 10 6 Pa. 
On the other hand, our simulations typically use a value of 
A ~ 0.01 LBU. This requires that the LBU of free energy 
density is about 10 8 Pa in SI units. Therefore our choice of 
elastic constant, K = 0.01 LBU, equates to K = 10" 10 N, 
corresponding to a Frank elastic constant K/2 = 50 pN which 
is sensible. 

For dynamical studies one can use the above results to cali- 
brate the lattice Boltzmann time unit, and then crosscheck that 
the resulting fluid density gives acceptable Reynolds numbers. 
This is indeed confirmed in [19] for our simulations, which 
choose an order parameter mobility F = 0.3 LBU in the no- 
tation of that paper, and a fluid viscosity r\ — 1 — 2 LBU. 
However, such considerations are not crucial in the current 
context where the dynamics need not be realistic, so long as 
they accurately minimize T. To summarize the results above 
(and the additional discussion in fT9ll ): our simulations faith- 
fully represent experimentally realisable BP-forming materi- 
als, subject to the interpretation of the LBU for length, energy 
density and time are close to 10 nm, 100 MPa and 1 ns respec- 
tively. 



FIG. 5: Initial (left) and final (right, after 300,000 lattice Boltzmann steps) configuration for a run which is started with randomly placed and 
oriented doubly twisted droplets, inside an isotropic background, as opposed to a cholesteric one as presented in the main text. Also in this 
case we observe a final amorphous network, very similar to the one shown in Fig. 1 of the main text. 



6. Free energy comparisons 

Table I shows the free energy obtained for our candidate 
blue phase III, blue phase II and O5 , in a temperature-chirality 
range somewhat larger than the one presented in the main 
text. (BPI is not competitive throughout this region.) At two 
points in parameter space, the free energies (in bold) are lower 
for BPIII than any other structure. For all other parameter 
sets in the Table, with the exception of one marked "CHOL", 
the BPIII candidate remains kinetically stable (when initiated 
from double-twist droplets in a cholesteric phase) even though 
it has higher free energy than one or more ordered phases. 

A difference in the tabulated values of Af = 0.1 be- 
tween phases (that is, Af = 10 -6 LBU) corresponds, with 
the parameter mapping detailed above, to a free energy dif- 
ference per unit cell of around per unit cell. This 
assumes a cell size of 160 nm, which is at the low end for 
blue phases; for a cell of say 500 nm the difference would be 
several hundreds ksT per cell. Adding thermal noise in the 
simulations (thereby weighting states by the Boltzmann factor 
exp[— F/ksT] rather than simply minimizing J 7 ) would be 
expected to create entropic shifts in the free energies of order 
fc#T per unit cell, reflecting configurational entropy and/or 
local order parameter fluctuations. These shifts would likely 



favor the amorphous BPIII structure, but, by these estimates, 
will do so only marginally. 

The above scaling of free energy densities via /c#T per unit 
cell is useful when discussing fluctuation effects but extrane- 
ous at the level of free energy minimization using the Landau 
- de Gennes functional, in which ksT does not appear. In that 
context, it is more natural to express the differences in free 
energy density between phases in units of either Aq or Kq$ . 
Since it relates to more easily measured quantities the latter 
scaling is preferable, and is followed in Fig. 2 of the main 
text. 



7. Dielectric anisotropy effect 

In the main text, we have considered the response of our 
candidate blue phase III structure to an electric field, assum- 
ing that the medium had a positive dielectric constant. We 
have also considered the case of negative dielectric constant, 
which favors an ordering of the director field normal to the 
electric field direction. Also in this case we observe ordering, 
although the final state is different, and consists of double he- 
lical disclinations parallel to the direction of the field. This 
state is shown in Fig. 6. 
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r 


K 


7 


A • 10 3 


r - BPIII 


f{r) - BPIII 


r-05 


f(r) - 05 


r - BPII 


f(r) - BPII 





2 


3 


6.9396 


1.097 


-3.104123 


1.009 


-3.047634 


0.908 


-3.161573 





2.5 


3 


4.4413 


1.090 


-1.804606 


1.016 


-1.828106 


0.914 


-1.791576 





2.75 


3 


3.6705 


1.049 


-1.415262 


1.018 


-1.465741 


0.917 


-1.376351 





3 


3 


3.0843 


1.097 


-1.144493 


1.020 


-1.192417 


0.919 


-1.062554 


-0.25 


1.75 


3.0857 


8.8122 


1.100 


-6.003057 


1.000 


-5.784407 


0.903 


-6.154037 


-0.25 


2 


3.0857 


6.7468 


1.103 


-4.42299 


1.007 


-4.292183 


0.907 


-4.508668 


-0.25 


2.25 


3.0857 


5.3308 


1.093 


-3.3654 


1.012 


-3.302209 


0.911 


-3.400327 


-0.25 


2.5 


3.0857 


4.3180 


1.085 


-2.632778 


1.015 


-2.597435 


0.914 


-2.602943 


-0.25 


2.6 


3.0857 


3.9922 


1.091 


-2.39750 


1.016 


-2.376423 


0.916 


-2.351800 


-0.25 


2.75 


3.0857 


3.5686 


1.043 


-2.08964 


1.017 


-2.092873 


0.917 


-2.046106 


-0.25 


3 


3.0857 


2.9986 


1.096 


-1.672261 


1.019 


-1.714777 


0.919 


-1.614218 


-0.5 


1.75 


3.1765 


8.5604 


1.026 


-5.101 (CHOL) 


0.996 


-7.662615 


0.901 


-8.166828 


-0.5 


2 


3.1765 


6.5540 


1.130 


-5.872775 


1.004 


-5.690254 


0.907 


-6.016180 


-0.5 


2.5 


3.1765 


4.1946 


1.069 


-3.501575 


1.013 


-3.472278 


0.913 


-3.546119 


-0.5 


2.75 


3.1765 


3.4666 


1.050 


-2.770071 


1.016 


-2.812765 


0.915 


-2.798638 


-0.5 


3 


3.1765 


2.9129 


1.045 


-2.250650 


1.018 


-2.319368 


0.918 


-2.2350 



TABLE I: Final free energy densities f(r) = 10 5 / in LBU. These are fully equilibrated (with redshift r as shown) for the ordered phases, 
and are redshifted also for the runs initialized from double twisted droplets in a cholesteric background for BPIII. Bold entries mark the free 
energies where BPIII is more stable than any of the ordered phases. (BPI is not competitive in this region of parameter space.) Note that for all 
parameters shown here BPIII is kinetically stable when initialized as described, with the exception of one low chirality run (marked CHOL) 
which reverted to a pure cholesteric state. 



